Multi-omics data elucidate parasite-host-microbiota interactions and resistance to Haemonchus contortus in sheep.

BACKGROUND: The integration of molecular data from hosts, parasites, and microbiota can enhance our understanding of the complex biological interactions underlying the resistance of hosts to parasites. Haemonchus contortus, the predominant sheep gastrointestinal parasite species in the tropics, causes significant production and economic losses, which are further compounded by the diminishing efficiency of chemical control owing to anthelmintic resistance. Knowledge of how the host responds to infection and how the parasite, in combination with microbiota, modulates host immunity can guide selection decisions to breed animals with improved parasite resistance. This understanding will help refine management practices and advance the development of new therapeutics for long-term helminth control. METHODS: Eggs per gram (EPG) of feces were obtained from Morada Nova sheep subjected to two artificial infections with H. contortus and used as a proxy to select animals with high resistance or susceptibility for transcriptome sequencing (RNA-seq) of the abomasum and 50 K single-nucleotide genotyping. Additionally, RNA-seq data for H. contortus were generated, and amplicon sequence variants (ASV) were obtained using polymerase chain reaction amplification and sequencing of bacterial and archaeal 16S ribosomal RNA genes from sheep feces and rumen content. RESULTS: The heritability estimate for EPG was 0.12. GAST, GNLY, IL13, MGRN1, FGF14, and RORC genes and transcripts were differentially expressed between resistant and susceptible animals. A genome-wide association study identified regions on chromosomes 2 and 11 that harbor candidate genes for resistance, immune response, body weight, and adaptation. Trans-expression quantitative trait loci were found between significant variants and differentially expressed transcripts. Functional co-expression modules based on sheep genes and ASVs correlated with resistance to H. contortus, showing enrichment in pathways of response to bacteria, immune and inflammatory responses, and hub features of the Christensenellaceae, Bacteroides, and Methanobrevibacter genera; Prevotellaceae family; and Verrucomicrobiota phylum. In H. contortus, some mitochondrial, collagen-, and cuticle-related genes were expressed only in parasites isolated from susceptible sheep. CONCLUSIONS: The present study identified chromosome regions, genes, transcripts, and pathways involved in the elaborate interactions between the sheep host, its gastrointestinal microbiota, and the H. contortus parasite. These findings will assist in the development of animal selection strategies for parasite resistance and interdisciplinary approaches to control H. contortus infection in sheep.

Expected values for the accuracy of predicted breeding values accounting for genetic differences between reference and target populations.

BACKGROUND: Genetic merit, or breeding values as referred to in livestock and crop breeding programs, is one of the keys to the successful selection of animals in commercial farming systems. The developments in statistical methods during the twentieth century and single nucleotide polymorphism (SNP) chip technologies in the twenty-first century have revolutionized agricultural production, by allowing highly accurate predictions of breeding values for selection candidates at a very early age. Nonetheless, for many breeding populations, realized accuracies of predicted breeding values (PBV) remain below the theoretical maximum, even when the reference population is sufficiently large, and SNPs included in the model are in sufficient linkage disequilibrium (LD) with the quantitative trait locus (QTL). This is particularly noticeable over generations, as we observe the so-called erosion of the effects of SNPs due to recombinations, accompanied by the erosion of the accuracy of prediction. While accurately quantifying the erosion at the individual SNP level is a difficult and unresolved task, quantifying the erosion of the accuracy of prediction is a more tractable problem. In this paper, we describe a method that uses the relationship between reference and target populations to calculate expected values for the accuracies of predicted breeding values for non-phenotyped individuals accounting for erosion. The accuracy of the expected values was evaluated through simulations, and a further evaluation was performed on real data. RESULTS: Using simulations, we empirically confirmed that our expected values for the accuracy of PBV accounting for erosion were able to correctly determine the prediction accuracy of breeding values for non-phenotyped individuals. When comparing the expected to the realized accuracies of PBV with real data, only one out of the four traits evaluated presented accuracies that were significantly higher than the expected, approaching h 2 . CONCLUSIONS: We defined an index of genetic correlation between reference and target populations, which summarizes the expected overall erosion due to differences in allele frequencies and LD patterns between populations. We used this correlation along with a trait's heritability to derive expected values for the accuracy ( R ) of PBV accounting for the erosion, and demonstrated that our derived E R | erosion is a reliable metric.

Feet and legs malformation in Nellore cattle: genetic analysis and prioritization of GWAS results.

Beef cattle affected by feet and legs malformations (FLM) cannot perform their productive and reproductive functions satisfactorily, resulting in significant economic losses. Accelerated weight gain in young animals due to increased fat deposition can lead to ligaments, tendon and joint strain and promote gene expression patterns that lead to changes in the normal architecture of the feet and legs. The possible correlated response in the FLM due to yearling weight (YW) selection suggest that this second trait could be used as an indirect selection criterion. Therefore, FLM breeding values and the genetic correlation between FLM and yearling weight (YW) were estimated for 295,031 Nellore animals by fitting a linear-threshold model in a Bayesian approach. A genome-wide association study was performed to identify genomic windows and positional candidate genes associated with FLM. The effects of single nucleotide polymorphisms (SNPs) on FLM phenotypes (affected or unaffected) were estimated using the weighted single-step genomic BLUP method, based on genotypes of 12,537 animals for 461,057 SNPs. Twelve non-overlapping windows of 20 adjacent SNPs explaining more than 1% of the additive genetic variance were selected for candidate gene annotation. Functional and gene prioritization analysis of candidate genes identified six genes (ATG7, EXT1, ITGA1, PPARD, SCUBE3, and SHOX) that may play a role in FLM expression due to their known role in skeletal muscle development, aberrant bone growth, lipid metabolism, intramuscular fat deposition and skeletogenesis. Identifying genes linked to foot and leg malformations enables selective breeding for healthier herds by reducing the occurrence of these conditions. Genetic markers can be used to develop tests that identify carriers of these mutations, assisting breeders in making informed breeding decisions to minimize the incidence of malformations in future generations, resulting in greater productivity and animal welfare.

deepGBLUP: joint deep learning networks and GBLUP framework for accurate genomic prediction of complex traits in Korean native cattle.

BACKGROUND: Genomic prediction has become widespread as a valuable tool to estimate genetic merit in animal and plant breeding. Here we develop a novel genomic prediction algorithm, called deepGBLUP, which integrates deep learning networks and a genomic best linear unbiased prediction (GBLUP) framework. The deep learning networks assign marker effects using locally-connected layers and subsequently use them to estimate an initial genomic value through fully-connected layers. The GBLUP framework estimates three genomic values (additive, dominance, and epistasis) by leveraging respective genetic relationship matrices. Finally, deepGBLUP predicts a final genomic value by summing all the estimated genomic values. RESULTS: We compared the proposed deepGBLUP with the conventional GBLUP and Bayesian methods. Extensive experiments demonstrate that the proposed deepGBLUP yields state-of-the-art performance on Korean native cattle data across diverse traits, marker densities, and training sizes. In addition, they show that the proposed deepGBLUP can outperform the previous methods on simulated data across various heritabilities and quantitative trait loci (QTL) effects. CONCLUSIONS: We introduced a novel genomic prediction algorithm, deepGBLUP, which successfully integrates deep learning networks and GBLUP framework. Through comprehensive evaluations on the Korean native cattle data and simulated data, deepGBLUP consistently achieved superior performance across various traits, marker densities, training sizes, heritabilities, and QTL effects. Therefore, deepGBLUP is an efficient method to estimate an accurate genomic value. The source code and manual for deepGBLUP are available at https://github.com/gywns6287/deepGBLUP .

The SWine IMputation (SWIM) haplotype reference panel enables nucleotide resolution genetic mapping in pigs.

Genetic mapping to identify genes and alleles associated with or causing economically important quantitative trait variation in livestock animals such as pigs is a major goal in animal genetic improvement. Despite recent advances in high-throughput genotyping technologies, the resolution of genetic mapping in pigs remains poor due in part to the low density of genotyped variant sites. In this study, we overcame this limitation by developing a reference haplotype panel for pigs based on 2259 whole genome-sequenced animals representing 44 pig breeds. We evaluated software combinations and breed composition to optimize the imputation procedure and achieved an average concordance rate in excess of 96%, a non-reference concordance rate of 88%, and an r2 of 0.85. We demonstrated in two case studies that genotype imputation using this resource can dramatically improve the resolution of genetic mapping. A public web server has been developed to allow the pig genetics community to fully utilize this resource. We expect this resource to facilitate genetic mapping and accelerate genetic improvement in pigs.

Runs of homozygosity and cross-generational inbreeding of Iranian fat-tailed sheep.

The Lori-Bakhtiari fat-tailed sheep is one of the most important heavyweight native breeds of Iran. The breed is robust and well-adapted to semi-arid regions and an important resource for smallholder farms. An established nucleus-based breeding scheme is used to improve their production traits but there is an indication of inbreeding depression and loss of genetic diversity due to selection. Here, we estimated the inbreeding levels and the distribution of runs of homozygosity (ROH) islands in 122 multi-generational female Lori-Bakhtiari from different half-sib families selected from a breeding station that were genotyped on the 50k array. A total of 2404 ROH islands were identified. On average, there were 19.70 ± 1.4 ROH per individual ranging between 6 and 41. The mean length of the ROH was 4.1 ± 0.14 Mb. There were 1999 short ROH of length 1-6 Mb and another 300 in the 6-12 Mb range. Additionally long ROH indicative of inbreeding were found in the ranges of 12-24 Mb (95) and 24-48 Mb (10). The average inbreeding coefficient (FROH) was 0.031 ± 0.003 with estimates varying from 0.006 to 0.083. Across generations, FROH increased from 0.019 ± 0.012 to 0.036 ± 0.007. Signatures of selection were identified on chromosomes 2, 6, and 10, encompassing 55 genes and 23 QTL associated with production traits. Inbreeding coefficients are currently within acceptable levels but across generations, inbreeding is increasing due to selection. The breeding program needs to actively monitor future inbreeding rates and ensure that the breed maintains or improves on its current levels of environmental adaptation.

Testicular hypoplasia in Nellore Cattle: Genetic analysis and functional analysis of genome-wide association study results.

Characterized by the incomplete development of the germinal epithelium of the seminiferous tubules, Testicular hypoplasia (TH) leads to decreased sperm concentration, increased morphological changes in sperm and azoospermia. Economic losses resulting from the disposal of affected bulls reduce the efficiency of meat production systems. A genome-wide association study and functional analysis were performed to identify genomic windows and the underlying positional candidate genes associated with TH in Nellore cattle. Phenotypic and pedigree data from 207,195 animals and genotypes (461,057 single nucleotide polymorphism, SNP) from 17,326 sires were used in this study. TH was evaluated as a binary trait measured at 18 months of age. A possible correlated response on TH resulting from the selection for scrotal circumference was evaluated by using a two-trait analysis. Thus, estimated breeding values were calculated by fitting a linear-threshold animal model in a Bayesian approach. The SNP effects were estimated using the weighted single-step genomic BLUP method. Twelve non-overlapping windows of 20 adjacent SNP that explained more than 1% of the additive genetic variance were selected for candidate gene annotation. Functional and gene prioritization analysis of the candidate genes identified genes (KHDRBS3, GPX5, STAR, ERLIN2), which might play an important role in the expression of TH due to their known roles in the spermatogenesis process, synthesis of steroids and lipid metabolism.

Evaluation of Whole-Genome Sequence Imputation Strategies in Korean Hanwoo Cattle.

This study evaluated the accuracy of sequence imputation in Hanwoo beef cattle using different reference panels: a large multi-breed reference with no Hanwoo (n = 6269), a much smaller Hanwoo purebred reference (n = 88), and both datasets combined (n = 6357). The target animals were 136 cattle both sequenced and genotyped with the Illumina BovineSNP50 v2 (50K). The average imputation accuracy measured by the Pearson correlation (R) was 0.695 with the multi-breed reference, 0.876 with the purebred Hanwoo, and 0.887 with the combined data; the average concordance rates (CR) were 88.16%, 94.49%, and 94.84%, respectively. The accuracy gains from adding a large multi-breed reference of 6269 samples to only 88 Hanwoo was marginal; however, the concordance rate for the heterozygotes decreased from 85% to 82%, and the concordance rate for fixed SNPs in Hanwoo also decreased from 99.98% to 98.73%. Although the multi-breed panel was large, it was not sufficiently representative of the breed for accurate imputation without the Hanwoo animals. Additionally, we evaluated the value of high-density 700K genotypes (n = 991) as an intermediary step in the imputation process. The imputation accuracy differences were negligible between a single-step imputation strategy from 50K directly to sequence and a two-step imputation approach (50K-700K-sequence). We also observed that imputed sequence data can be used as a reference panel for imputation (mean R = 0.9650, mean CR = 98.35%). Finally, we identified 31 poorly imputed genomic regions in the Hanwoo genome and demonstrated that imputation accuracies were particularly lower at the chromosomal ends.

Genome Wide Scan to Identify Potential Genomic Regions Associated With Milk Protein and Minerals in Vrindavani Cattle.

In this study, genome-wide association study (GWAS) was conducted for identifying significantly associated genomic regions/SNPs with milk protein and minerals in the 96 taurine-indicine crossbred (Vrindavani) cows using 50K SNP Chip. After quality control, a total of 41,427 SNPs were retained and were further analyzed using a single-SNP additive linear model. Lactation stage, parity, test day milk yield and proportion of exotic inheritance were included as fixed effects in GWAS model. Across all traits, 13 genome-wide significant (p < 1.20 x 10-06) and 49 suggestive significant (p < 2.41 x 10-05) SNPs were identified which were located on 18 different autosomes. The strongest association for protein percentage, calcium (Ca), phosphorus (P), copper (Cu), zinc (Zn), and iron (Fe) were found on BTA 18, 7, 2, 3, 14, and 2, respectively. No significant SNP was detected for manganese (Mn). Several significant SNPs identified were within or close proximity to CDH13, BHLHE40, EDIL3, HAPLN1, INHBB, USP24, ZFAT, and IKZF2 gene, respectively. Enrichment analysis of the identified candidate genes elucidated biological processes, cellular components, and molecular functions involved in metal ion binding, ion transportation, transmembrane protein, and signaling pathways. This study provided a groundwork to characterize the molecular mechanism for the phenotypic variation in milk protein percentage and minerals in crossbred cattle. Further work is required on a larger sample size with fine mapping of identified QTL to validate potential candidate regions.

Long-term artificial selection of Hanwoo (Korean) cattle left genetic signatures for the breeding traits and has altered the genomic structure.

Indigenous Korean breeds such as Hanwoo (Korean) cattle have adapted to their local environment during the past 5000 years. In the 1980s, the National Genetic Improvement Program was established to develop a modern economic breed for beef production in Korea through artificial selection. This process is thought to have altered the genomic structure of breeding traits over time. The detection of genetic variants under selection could help to elucidate the genetic mechanism of artificial selection in modern cattle breeds. Indigenous Hanwoo cattle have adapted in response to local natural and artificial selection during a 40-year breeding program. We analyzed genomic changes in the selection signatures of an unselected population (USP; n = 362) and a selected population (KPN; n = 667) of Hanwoo cattle. Genomic changes due to long-term artificial selection were identified using a genome-wide integrated haplotype score (iHS) and a genome-wide association study (GWAS). Signatures of recent selection were detected as positive (piHS > 6) or negative (piHS < -6) iHS scores spanning more than 46 related genes in KPN cattle, but none in USP cattle. A region adjacent to the PLAG1 gene was found to be under strong selection for carcass weight. The GWAS results also showed a selection signature on BTA14, but none on BTA13. Pathway and quantitative trait locus analysis results identified candidate genes related to energy metabolism, feed efficiency, and reproductive traits in Hanwoo cattle. Strong selection significantly altered Hanwoo cattle genome structural properties such as linkage disequilibrium (LD) and haplotypes through causal mutation for target traits. Haplotype changes of genome structure which are changes of ancestral allele to derived alleles due to selection were clearly identified on BTA13 and BTA14; however, the structure of the LD block was not clearly observed except BTA14. Thus, selection based on EBVs would be working very well in Hanwoo cattle breeding program appears to have been highly successful.

GPS Coordinates for Modelling Correlated Herd Effects in Genomic Prediction Models Applied to Hanwoo Beef Cattle.

It is widely known that the environment influences phenotypic expression and that its effects must be accounted for in genetic evaluation programs. The most used method to account for environmental effects is to add herd and contemporary group to the model. Although generally informative, the herd effect treats different farms as independent units. However, if two farms are located physically close to each other, they potentially share correlated environmental factors. We introduce a method to model herd effects that uses the physical distances between farms based on the Global Positioning System (GPS) coordinates as a proxy for the correlation matrix of these effects that aims to account for similarities and differences between farms due to environmental factors. A population of Hanwoo Korean cattle was used to evaluate the impact of modelling herd effects as correlated, in comparison to assuming the farms as completely independent units, on the variance components and genomic prediction. The main result was an increase in the reliabilities of the predicted genomic breeding values compared to reliabilities obtained with traditional models (across four traits evaluated, reliabilities of prediction presented increases that ranged from 0.05 ± 0.01 to 0.33 ± 0.03), suggesting that these models may overestimate heritabilities. Although little to no significant gain was obtained in phenotypic prediction, the increased reliability of the predicted genomic breeding values is of practical relevance for genetic evaluation programs.

Genomic Prediction for Whole Weight, Body Shape, Meat Yield, and Color Traits in the Portuguese Oyster Crassostrea angulata.

Genetic improvement for quality traits, especially color and meat yield, has been limited in aquaculture because the assessment of these traits requires that the animals be slaughtered first. Genotyping technologies do, however, provide an opportunity to improve the selection efficiency for these traits. The main purpose of this study is to assess the potential for using genomic information to improve meat yield (soft tissue weight and condition index), body shape (cup and fan ratios), color (shell and mantle), and whole weight traits at harvest in the Portuguese oyster, Crassostrea angulata. The study consisted of 647 oysters: 188 oysters from 57 full-sib families from the first generation and 459 oysters from 33 full-sib families from the second generation. The number per family ranged from two to eight oysters for the first and 12-15 oysters for the second generation. After quality control, a set of 13,048 markers were analyzed to estimate the genetic parameters (heritability and genetic correlation) and predictive accuracy of the genomic selection for these traits. The multi-locus mixed model analysis indicated high estimates of heritability for meat yield traits: 0.43 for soft tissue weight and 0.77 for condition index. The estimated genomic heritabilities were 0.45 for whole weight, 0.24 for cup ratio, and 0.33 for fan ratio and ranged from 0.14 to 0.54 for color traits. The genetic correlations among whole weight, meat yield, and body shape traits were favorably positive, suggesting that the selection for whole weight would have beneficial effects on meat yield and body shape traits. Of paramount importance is the fact that the genomic prediction showed moderate to high accuracy for the traits studied (0.38-0.92). Therefore, there are good prospects to improve whole weight, meat yield, body shape, and color traits using genomic information. A multi-trait selection program using the genomic information can boost the genetic gain and minimize inbreeding in the long-term for this population.

Meta-Analysis of Heifer Traits Identified Reproductive Pathways in Bos indicus Cattle.

Fertility traits measured early in life define the reproductive potential of heifers. Knowledge of genetics and biology can help devise genomic selection methods to improve heifer fertility. In this study, we used ~2400 Brahman cattle to perform GWAS and multi-trait meta-analysis to determine genomic regions associated with heifer fertility. Heifer traits measured were pregnancy at first mating opportunity (PREG1, a binary trait), first conception score (FCS, score 1 to 3) and rebreeding score (REB, score 1 to 3.5). The heritability estimates were 0.17 (0.03) for PREG1, 0.11 (0.05) for FCS and 0.28 (0.05) for REB. The three traits were highly genetically correlated (0.75-0.83) as expected. Meta-analysis was performed using SNP effects estimated for each of the three traits, adjusted for standard error. We identified 1359 significant SNPs (p-value < 9.9 × 10-6 at FDR < 0.0001) in the multi-trait meta-analysis. Genomic regions of 0.5 Mb around each significant SNP from the meta-analysis were annotated to create a list of 2560 positional candidate genes. The most significant SNP was in the vicinity of a genomic region on chromosome 8, encompassing the genes SLC44A1, FSD1L, FKTN, TAL2 and TMEM38B. The genomic region in humans that contains homologs of these genes is associated with age at puberty in girls. Top significant SNPs pointed to additional fertility-related genes, again within a 0.5 Mb region, including ESR2, ITPR1, GNG2, RGS9BP, ANKRD27, TDRD12, GRM1, MTHFD1, PTGDR and NTNG1. Functional pathway enrichment analysis resulted in many positional candidate genes relating to known fertility pathways, including GnRH signaling, estrogen signaling, progesterone mediated oocyte maturation, cAMP signaling, calcium signaling, glutamatergic signaling, focal adhesion, PI3K-AKT signaling and ovarian steroidogenesis pathway. The comparison of results from this study with previous transcriptomics and proteomics studies on puberty of the same cattle breed (Brahman) but in a different population identified 392 genes in common from which some genes-BRAF, GABRA2, GABR1B, GAD1, FSHR, CNGA3, PDE10A, SNAP25, ESR2, GRIA2, ORAI1, EGFR, CHRNA5, VDAC2, ACVR2B, ORAI3, CYP11A1, GRIN2A, ATP2B3, CAMK2A, PLA2G, CAMK2D and MAPK3-are also part of the above-mentioned pathways. The biological functions of the positional candidate genes and their annotation to known pathways allowed integrating the results into a bigger picture of molecular mechanisms related to puberty in the hypothalamus-pituitary-ovarian axis. A reasonable number of genes, common between previous puberty studies and this study on early reproductive traits, corroborates the proposed molecular mechanisms. This study identified the polymorphism associated with early reproductive traits, and candidate genes that provided a visualization of the proposed mechanisms, coordinating the hypothalamic, pituitary, and ovarian functions for reproductive performance in Brahman cattle.

Identification of genes affecting milk fat and fatty acid composition in Vrindavani crossbred cattle using 50 K SNP-Chip.

The aim of this study was to identify candidate genes associated with milk fat per cent and fatty acid (FA) composition in Vrindavani cattle using the Illumina 50 K single-nucleotide polymorphism (SNP) array. After quality control, a total of 41,427 informative and high-quality SNPs were used for a genome-wide association study (GWAS) for milk fat percentage and 16 different types of fatty acids. Lactation stage, parity, test day milk yield, and proportion of exotic inheritance were included as fixed effects in the GWAS model. A total of 67 genome-wide significant (P < 1.20 × 10-06) SNPs and 176 suggestive significant (P < 2.41 × 10-05) SNPs were identified. Out of these, 15 SNPs were associated with more than one trait. The strongest associations were found on BTA14 for milk fat percentage and on BTA2 and BTA16 for polyunsaturated fatty acids. Several significant SNPs were identified close to or within the genes ELOVL6, FABP4, PMP2, PLIN1, MFGE8, GHRL2, and LDLRAD3 which are known to be associated with fat percentage and FA composition in dairy cattle breeds. This study is a step forward to better characterize the molecular mechanisms of phenotypic variation in milk fatty acids in a taurine-indicine composite cattle breed reared in tropical environments.

Heuristic hyperparameter optimization of deep learning models for genomic prediction.

There is a growing interest among quantitative geneticists and animal breeders in the use of deep learning (DL) for genomic prediction. However, the performance of DL is affected by hyperparameters that are typically manually set by users. These hyperparameters do not simply specify the architecture of the model; they are also critical for the efficacy of the optimization and model-fitting process. To date, most DL approaches used for genomic prediction have concentrated on identifying suitable hyperparameters by exploring discrete options from a subset of the hyperparameter space. Enlarging the hyperparameter optimization search space with continuous hyperparameters is a daunting combinatorial problem. To deal with this problem, we propose using differential evolution (DE) to perform an efficient search of arbitrarily complex hyperparameter spaces in DL models, and we apply this to the specific case of genomic prediction of livestock phenotypes. This approach was evaluated on two pig and cattle datasets with real genotypes and simulated phenotypes (N = 7,539 animals and M = 48,541 markers) and one real dataset (N = 910 individuals and M = 28,916 markers). Hyperparameters were evaluated using cross-validation. We compared the predictive performance of DL models using hyperparameters optimized by DE against DL models with "best practice" hyperparameters selected from published studies and baseline DL models with randomly specified hyperparameters. Optimized models using DE showed a clear improvement in predictive performance across all three datasets. DE optimized hyperparameters also resulted in DL models with less overfitting and less variation in predictive performance over repeated retraining compared to non-optimized DL models.

Prediction Accuracies of Genomic Selection for Nine Commercially Important Traits in the Portuguese Oyster (Crassostrea angulata) Using DArT-Seq Technology.

Genomic selection has been widely used in terrestrial animals but has had limited application in aquaculture due to relatively high genotyping costs. Genomic information has an important role in improving the prediction accuracy of breeding values, especially for traits that are difficult or expensive to measure. The purposes of this study were to (i) further evaluate the use of genomic information to improve prediction accuracies of breeding values from, (ii) compare different prediction methods (BayesA, BayesCπ and GBLUP) on prediction accuracies in our field data, and (iii) investigate the effects of different SNP marker densities on prediction accuracies of traits in the Portuguese oyster (Crassostrea angulata). The traits studied are all of economic importance and included morphometric traits (shell length, shell width, shell depth, shell weight), edibility traits (tenderness, taste, moisture content), and disease traits (Polydora sp. and Marteilioides chungmuensis). A total of 18,849 single nucleotide polymorphisms were obtained from genotyping by sequencing and used to estimate genetic parameters (heritability and genetic correlation) and the prediction accuracy of genomic selection for these traits. Multi-locus mixed model analysis indicated high estimates of heritability for edibility traits; 0.44 for moisture content, 0.59 for taste, and 0.72 for tenderness. The morphometric traits, shell length, shell width, shell depth and shell weight had estimated genomic heritabilities ranging from 0.28 to 0.55. The genomic heritabilities were relatively low for the disease related traits: Polydora sp. prevalence (0.11) and M. chungmuensis (0.10). Genomic correlations between whole weight and other morphometric traits were from moderate to high and positive (0.58-0.90). However, unfavourably positive genomic correlations were observed between whole weight and the disease traits (0.35-0.37). The genomic best linear unbiased prediction method (GBLUP) showed slightly higher accuracy for the traits studied (0.240-0.794) compared with both BayesA and BayesCπ methods but these differences were not significant. In addition, there is a large potential for using low-density SNP markers for genomic selection in this population at a number of 3000 SNPs. Therefore, there is the prospect to improve morphometric, edibility and disease related traits using genomic information in this species.

Differential Gene Expression in Longissimus Dorsi Muscle of Hanwoo Steers-New Insight in Genes Involved in Marbling Development at Younger Ages.

The Korean Hanwoo breed possesses a high capacity to accumulate intramuscular fat, which is measured as a marbling score in the beef industry. Unfortunately, the development of marbling is not completely understood and the identification of differentially expressed genes at an early age is required to better understand this trait. In this study, we took muscle samples from 12 Hanwoo steers at the age of 18 and 30 months. From the contrast between age and marbling score, we identified in total 1883 differentially expressed genes (FDR < 0.05 and logarithm fold change ≥ 1.5) with 782 genes up-regulated and 1101 down-regulated. Differences in gene expression were higher between the ages x marbling groups rather than between high and low marbling groups. At 18 months of age, the genes SLC38A4, ABCA10, APOL6, and two novel genes (ENSBTAG00000015330 and ENSBTAG00000046041) were up-regulated in the high marbling group. From the protein-protein interaction network analysis, we identified unique networks when comparing marbling scores between different ages. Nineteen genes (AGT, SERPINE1, ADORA1, FOS, LEP, FOXO1, FOXO3, ADIPOQ, ITGA1, SDC1, SDC4, ITGB3, ITGB4, CXCL10, ACTG2, MX1, EDN1, ACTA2, and ESPL1) were identified to have an important role in marbling development. Further analyses are needed to better understand the role of these genes.

Signatures of Selection in Composite Vrindavani Cattle of India.

Vrindavani is an Indian composite cattle breed developed by crossbreeding taurine dairy breeds with native indicine cattle. The constituent breeds were selected for higher milk production and adaptation to the tropical climate. However, the selection response for production and adaptation traits in the Vrindavani genome is not explored. In this study, we provide the first overview of the selection signatures in the Vrindavani genome. A total of 96 Vrindavani cattle were genotyped using the BovineSNP50 BeadChip and the SNP genotype data of its constituent breeds were collected from a public database. Within-breed selection signatures in Vrindavani were investigated using the integrated haplotype score (iHS). The Vrindavani breed was also compared to each of its parental breeds to discover between-population signatures of selection using two approaches, cross-population extended haplotype homozygosity (XP-EHH) and fixation index (F ST). We identified 11 common regions detected by more than one method harboring genes such as LRP1B, TNNI3K, APOB, CACNA2D1, FAM110B, and SPATA17 associated with production and adaptation. Overall, our results suggested stronger selective pressure on regions responsible for adaptation compared to milk yield.

Prediction of marbling score and carcass traits in Korean Hanwoo beef cattle using machine learning methods and synthetic minority oversampling technique.

Pricing of Hanwoo beef in the Korean market is primarily based on meat quality, and particularly on marbling score. The ability to accurately predict marbling score early in the life of an animal is extremely valuable for producers to meet the requirements of their target market, and for genetic selection. A total of 3989 Korean Hanwoo cattle (2108 with 50 k SNP genotypes) and 45 phenotypic features were available for this study. Four machine learning (ML) algorithms were applied to predict six carcass traits and compared against linear regression prediction models. In most scenarios, SMO was the best performing algorithm. The most and least accurately predicted traits were carcass weight and marbling score with correlation of 0.95 and 0.64 respectively. Additionally, the value of using a synthetic minority over-sampling technique (SMOTE) was evaluated and results showed a slight improvement in the prediction error of marbling score. Machine Learning approaches can be useful tools to predict important carcass traits in beef cattle.